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ABSTRACT 

The detection of bright X-ray flares superimposed on the regular afterglow decay 
in Swift gamma-ray bursts has triggered theoretical speculations on their origin. We 
study the temporal properties of flares due to internal dissipation and external shock 
mechanisms. We first show that at least a sizable fraction of the flares cannot be related 
to external shock mechanisms, since external shock flares evolve on much longer time 
scales than observed. We then study flares from internal dissipation, showing that the 
temporal properties allow us to distinguish the emission of slow early shells from that 
of late faster shells. We show that, due to the rapid evolution of the detected flares, 
it is most likely that the flares are produced by relatively fast shells ejected by the 
central engine shortly before they are observed. This implies that the central engine 
must be active for, in some cases, as long as one day. We finally discuss the constraints 
and implications that this observation has on the properties and physics of the inner 
engine, and we elaborate on possible future observational tests on the flare sample to 
further understand their origin and physics. 
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1 INTRODUCTION 

The external shock model predicts afterglows characterized 
by a power-law decay in time, with the possible presence of 
breaks connecting branches of different slope (Meszaros & 
Rees 1997). Such changes are due to either geometrical prop- 
erties of the fireball (geometrical beaming, Rhoads 1999) or 
to spectral transitions (e.g. from a cooling to a non-cooling 
electron population. Sari, Piran & Narayan 1998). In the 
pre- Swift era, most afterglows were consistent with the sim- 
plest version of this model. There were, however, notable 
exceptions, such as GRB 000301C (Masetti et al. 2000), 
GRB 021004 (Lazzati et al. 2002), and GRB 030329 (Mathe- 
son et al. 2003). All these bursts displayed optical variability 
in the form of bumps or wiggles superimposed on the smooth 
power-law decay. 

A number of explanations have been discussed in the lit- 
erature to account for the variations in the optical afterglow 
brightness. These include inhomogeneities in the external 
density (Wang & Loeb 2000; Lazzati et al. 2002; Heyl & 
Perna 2003), refreshed shocks due to the collision of a late 
shell of plasma with the external shock material (Rees & 
Meszaros 1998), angular inhomogeneities in the fireball en- 
ergy distribution (Nakar, Piran & Granot 2002) and grav- 
itational lensing (Loeb & Perna 1998; Garnavich, Loeb & 
Stanek 2000). Different mechanisms to produce variability 
in afterglows can in principle be distinguished through their 



temporal and spectral properties. In practice, however, a 
consensus has not been yet reached due to the lack of un- 
ambiguous observations. 

More recently. Swift observations have revealed that 
fiaring activity is relatively common in the early phases of 
GRB afterglows, sometimes extending for over a day (e.g., 
Chincarini 2006; O'Brien et al. 2006). The temporal prop- 
erties of the fiares, their intensity, and their spectra suggest 
an origin that is unrelated to the external shock (Chincar- 
ini 2006; Falcone et al. 2006), at least for a fraction of the 
bumps. Late episodes of "prompt emission'j^ can have two 
different origin. One possibility is that the inner engine it- 
self is active for a time as long as the detection time of the 
X-ray fiare. Alternatively, the engine can be short-lived but 
produce, together with the fast ejecta, a tail of slower ma- 
terial. Such slower material can produce internal dissipation 
(and therefore prompt emission) at late time. Slow shells 
ejected immediately after the fast ejecta will not produce 
an external shock at the canonical external shock radius 
since the ambient medium has been swept by the external 
shock produced by the fast ejecta, which develops earlier. 



^ Here and in the following we will call prompt emission all the 
radiation produced by the outflow inside the external shock. Re- 
lated mechanisms include, but are not limited to, internal shocks, 
magnetic dissipation, and comptonization of external photons. 
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Figure 1. Cartoon showing the geometry of the angular time 
scale. 



Distinguishing between these two scenarios bears important 
imphcations for the physics of the GRB engine. 

In this paper we discuss the timescales of flares in the 
different scenarios. We show that at least a large fraction of 
X-ray flares are due to a long lasting activity of the GRB 
engine, rather than to external shock activity or to the emis- 
sion of slow shells immediately after the prompt emission 
ends. 

This paper is organized as follows: in Sect. 2 we com- 
pute external shock timescales while in Sect. 3 we consider 
prompt emission flares and in Sect. 4 we discuss our results 
and compare them to Swift data. We summarize our flndings 
in Sect. 5. 



2 EXTERNAL SHOCK 

Consider flaring activity produced by a sudden brightening 
of the external shock. We consider now a brightening that 
involves the whole surface of the external shock. Even if the 
duration of the activity in the comoving frame is negligible, 
the observer at inflnity still detects photons over a finite 
amount of time, due to the "curvature effect" (see Fig. [T|. 
The shape of the observed pulse is therefore the narrowest 
possible in time. A longer activity in the comoving frame 
will produce a longer pulse, obtained by convolution of the 
two functional shapes. 

Let us define the quantity E'{v') to be the energy re- 
leased by the unit area of the fireball as a result of the flaring 
activity, in the comoving frame. The flux received at infinity 
is given by (see also Kumar & Panaitescu 2000): 



E'{v/5) 52 dS 



where D is the distance to the source, 5 ; 
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Figure 2. Minimum time scale ratio for flares in different dy- 
namical settings as a function of the spectral index rj. From top 
to bottom, we show a radiative shock propagating into a uniform 
medium (a = 3), an adiabatic shock propagating into a uniform 
medium {a = 3/2), a radiative shock propagating into a wind 
profile (0 = 1), an adiabatic shock propagating into a wind pro- 
file {a = 1/2), and the limit for internal dissipation from a slow 
shell ejected during or immediately after the prompt phase. The 
shaded area represents the observational distribution oi At/t from 
Chincarini (2006), corrected for our At definition. 



is the Doppler factor, E is the fireball surface and tobs the 
time in the observer frame. This is given by: 
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Equation ([2} holds for any self-similar dynamical evo- 
lution where Fgh oc as long as Fsh ^ 1. Note that we 
indicate with the subscript -sii" quantities (such as Lorentz 
factor and speed) of the external shock, and without sub- 
script the same quantities for the material just behind the 
shock. It can be shown that Fgh ~ \/2F (e.g., Sari 1997). 

Easy considerations show that the only non constant 
quantity in Eq.([T} is the Doppler factor, since dE/dtobs = 
27r R c, and all the radiation is emitted at constant radius R. 
If we assume that we are far from a break in the afterglow 
spectrum, E'{v') oc and the functional shape of Fv{t) 

can be rewritten as: 



Solving Eq. ^ for cos 9, 
the start time of the flare to 
proflle can be written as: 



F^t) oc 



-('7+2) 



(3) 

and expressing time in unit of 
= i?/[2cF^i,(2a + 1)], the flare 



(4) 



where r = t/tg. 

Equation|4]can be used to provide a lower limit to an ob- 
servationally sound deflnition of the flare duration. Deflning 
the duration as the Full Width at Half Maximum (FWHM) 
of the pulse, and the reference time as the time tpk of the 
flare maximum, we obtain: 
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T 



> ( 2^ - 1 (4a + 2) 



(5) 



Equation ([5]) depends on both the dynamics and spec- 
trum of the external shock. A steeper spectrum and/or a 
more slowly evolving fireball will cause a narrower flare. 
Note that in this case (i.e. impulsive rebrightening of the 
shock), the maximum coincides with the start time. Should 
the activity in the comoving frame last for a finite amount 
of time, a rising phase would be observed. A more detailed 
treatment, including the width of the shell of shocked ma- 
terial would also provide a profile with a finite rising time. 
Figure [2] shows the time scale ratio for four typical values of 
Q as a function of the spectral index. The most relevant case 
for X-ray flares is the adiabatic ISM case, with a spectral 
slope rj ~ 1.2. We find for that specific case: At/t ^ 2. 

In an analogous context, Zhang et al. (2006) discussed 
the possible production of flares by various mechanisms such 
as density bumps, post-energy injection in the blast wave, 
and two-components or patchy jets. In all cases, they con- 
cluded that the characteristics of the flares were not gen- 
erally consistent with what predicted by these models (i.e. 
At/t ~ 1 and a slower-than observed flux decline). Here, 
we have derived detailed, quantitative, upper limits on the 
value of A/t. Our constraints are somewhat tighter than 
what previously derived by different authors in the context 
of the external model for flares (e.g., loka, Kobayashi & 
Zhang 2005) but comparable to more recent results (Wu 
et al. 2006). Our method differs in two important ways from 
what done in previous work. First, we have computed the ac- 
tual functional shape of the flare, and defined rigorously the 
time reference we adopt (the peak) and the time width (the 
FWHM). In previous works, these quantities were usually 
poorly defined, and the results could only be approximate. 
Secondly, we considered the difference between the Lorentz 
factor of the shock and that of the shocked material. This 
factor by itself leads to flares which are broader by a factor 
of two. 

Chincarini (2006) presents a distribution of At/t for a 
sample of Smft X-ray flares that is defined as the ratio of the 
Gaussian a over the peak time. The distribution is shown 
in their Fig. 8 and shaded in Fig. [5] Since for a Gaussian 
the FWHM is equal to 2.35(t, we can easily convert their 
definition of At into our definition. We conclude that the 
distribution of At/t of Swtft X-ray flares is characterized by 
0.02 ^ At/t ^ 1.3 and has average At/t — 0.3. Considering 
Fig. [2j we conclude that most, if not all, the detected re- 
brightenings are not related to events taking place on the 
external shock, unless only a small portion of the shock is 
involved in the rebrightening. Flares caused by a small por- 
tion of the shock are possible, but are unlikely to be very 
strong. Nakar & Granot (2006) studied the effect of a blob 
of high density interacting with the shock. They find that 
the rebrightening is minor, and characterized by a very slow 
rise. Alternatively, a large rebrightening on a small timescale 
can be observed if a narrow opening angle shell refreshes a 
small part of the external shock (Granot, Nakar & Piran 
2003). This is not likely to be the case for the Swift X-ray 
flares, since most of them are observed before the jet break 
time (Chincarini 2006, O'Brien 2006) and the opening an- 
gle of the jet grows with time (Lazzati & Begelman 2005; 
Morsony, Lazzati & Begelman 2006) rather than decrease. 



3 PROMPT EMISSION 

Since the largest majority of flares is unlikely to be pro- 
duced within the external shock, it must be produced dur- 
ing the initial phase of internal dissipation. In the following 
we therefore discuss flaring activity in the context of the 
prompt emission. In this case one can envisage two possi- 
ble scenarios for their production: (a) shells that are pro- 
duced during the GRB phase but that dissipate at much 
later times; (b) shells produced at late times by a long-lived 
engine and dissipating on a timescale comparable to that 
of the engine duration. Zhang et al. (2006) also concluded 
that the flares have an internal origin, but directly assumed 
that the engine must have been long-lived. In this section we 
try to discriminate between the scenarios (a) and (b). Being 
able to discriminate between these two scenarios bears im- 
portant implications for our understanding of the physics of 
the inner GRB engine (see more discussion on this in §4). 
In the following, we discuss the two scenarios above, pro- 
vide diagnostics for each of them, and show that the current 
Swift data is already able to discriminate between the two 
flare-production mechanisms. 

(a) Flaring activity due to dissipation within a freely 
expanding flow, released during the prompt GRB phase or 
shortly afterwards. This case has two important differences 
with respect to the external shock case studied above. First, 
there is no propagating shock, and therefore there is no dis- 
tinction between F and Fsh. Second, the flare is produced 
by material that has been coasting at constant F, rather 
than by material that has been slowing down. These two 
differences change considerably the equations above. 

Without repeating the derivation of Sect. 2, it can be 
easily seen that the result for the freely expanding flow can 
be obtained by substituting 4a -|- 2 with 2a -1-1 and adopting 
a = in Eq. [1 This yields: 



2~ 



l-pk 



(6) 



which, for the typical case of ~ 1.2 yields At/t ^ 0.25 (see 
Wu et al. 2006 for a similar computation restricted to the 
case of internal shocks). This number is tantalizing similar 
to the average width of Swift X-ray flares (Chincarini 2006; 
Fig. [21 . Yet, there is a sizable number of events for which 
the condition is violated by almost an order of magnitude. In 
addition, should the condition be barely satisfied, the decay 
part of the flare would be dominated by large angle emission 
(Kumar & Panaitescu 2000; Liang et al. 2006). This seems 
not to be the case, at least in some of the flares (Guetta et 
al. 2006). Note however that the large angle emission should 
not dominate any flare with a duration significantly larger 
than the one given in Eq.([6ll. 

(b) Shell that is ejected from the central engine after a 
time taj sizably larger than the prompt emission timescale 
Tgo- If radiation is released at radius R, it will be observed 
at a time (cfr. Eq. [2]) 



tobs = tcj -f — I 1 - cos 6^ 4- 



(7) 



which produces a shortening of the observed variability 
timescale 



At 



(8) 
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Equation (|8} allows for any arbitrarily small time scale, 
as long as the shell that produces radiation is ejected at a 
time comparable to the detection time. Therefore, the detec- 
tion of fast flares implies that the inner engine is active for 
a time much longer than the canonical prompt phase (Tgo). 

The observation of a flare that peaks at time tpk and is 
characterized by a width At allows us to constrain the time 
at which it was ejected from the central engine: 



-L ^ — ~ — ^ -L . yu, 

2 2+Tj \ i-pk 

4 DISCUSSION 

Equations (|6} and ([8| are general and hold for any 
dissipation and radiation mechanism that take place in 
the whole jet. Locally beamed emission mechanisms can 
produce fast variability (see, e.g., Lyutikov & Bland- 
ford 2004), but the peak luminosity of the flares would 
drop faster than observed as a function of tpk. Con- 
sider now the flares observed by Swift (for a recent 
review, see http://www.swift.ac.uk/rs06/Burrows.pdf). 
They can be observed up to ~ lO'^ s after the GRB onset 
and are characterized by a FWHM time 5t ~ 0.3 1. Applying 
this observational constraint implies that at least in about 
half of the cases the ejection time of the material from the 
central engine is comparable to the time at which the flare 
is observed. In other words, GRB engines are active well be- 
yond the observed Tgo of the prompt emission, in some cases 
for up to ~ 10^ s or a day timescale. 

This provides an important constraint to the properties 
of the inner engine of long duration GRBs. More can be 
obtained by further analysis of the properties of the flares. 
Flares are observed over a wide range of times, spanning 
about three orders of magnitude (from 100 to 10^ s after 
the end of the prompt emission). The distribution of 5t/t is 
however much narrower, spanning approximately only one 
order of magnitude. This implies a correlation between the 
ejection time and the duration of the ejection episode. The 
longer is the time after the engine turns on, the longer is 
the ejection episode. Interestingly, an analogous property 
was noted in GRBs with multiple prompt emission episodes 
(Ramirez- Ruiz & Merloni 2001). During the prompt phase, 
the duration of an emission episode is strongly correlated 
with the length of the quiescent time that precedes it. In 
the case of flares we can conclude that the correlation is 
with the total time since the engine onset; however, this 
does not imply the absence of a correlation with the length 
of the quiescent time. 

Our results bear important implications for an under- 
standing of the physical processes governing the GRB in- 
ner engine, which provides the ultimate source of energy to 
power GRBs. The GRB "inner engine" is believed to be a hy- 
peraccreting accretion disk. In the coUapsar model (Paczyn- 
sky 1998; MacFadyen & Woosley 1999), this is formed by the 
fallback material from the collapsing envelope of the mas- 
sive star, and the timescale for the duration of the accretion 
episode is set by the dynamical timescale of the collapsing 
envelope itself, which is on the order of several tens of sec- 
onds. On the other hand, in the binary merger model (e.g. 
neutron star-neutron star or neutron star-black hole; Eichler 
et al. 1989) , the accretion material is provided by the debris 



of the tidally disrupted neutron star (or white dwarf). In 
this case, the timescale over which the available material is 
accreted is set by the viscous timescale of the disk, which is 
a fraction of a second. 

The accretion timescale in the coUapsar model easily 
accounts for the duration of the prompt emission for the 
class of long bursts, while the accretion timescale in the 
binary merger model naturally yields the required timescales 
for the prompt emission of short bursts. Several types of 
observations in the last few years have indeed provided a 
strong support for the association between long GRBs and 
the collapse of a massive star (Stanek et al. 2003; Hjorth 
et al. 2003), while several pieces of evidence are gradually 
mounting toward the association of compact-object mergers 
and short bursts (Fox et al. 2005; Bloom et al. 2006; Nakar, 
Gal- Yam & Fox 2006). 

In the scenarios above, a relativistic outflow is ejected 
from the central engine, and radiation is produced as the 
bulk kinetic energy is dissipated and radiated at a given ra- 
dius. There is no a-priori constraint on the distribution of 
the Lorentz factor in the outflow. If the Lorentz factor of 
the ejecta should decrease with the ejection time, late time 
internal activity could be seen without the need for a long- 
living engine. We have shown that this is not the case, at 
least in a sizable fraction of the flares. This result makes the 
studies of possible mechanisms of reactivation of the GRB 
engine especially timely. King et al. (2005) suggested that 
fragmentation of the collapsing star could explain a (sin- 
gle) flare in the case of long bursts, while Dai et al. (2006) 
proposed a new mechanism of binary merger that can in- 
corporate the presence of a flare for short bursts. Perna, 
Armitage & Zhang (2006) noted that the fact that flares 
are observed in both long and short bursts, and with similar 
characteristics in both cases, is suggestive that the place of 
formation be in what is common between the two classes 
of bursts, i.e. the accretion disk. A disk that fragments as 
a result of gravitational instabilities in its outer parts, as 
supported by studies of hyperaccreting disks (Di Matteo, 
Perna & Narayan 2002; Chen & Beloborodov 2006), can 
account for flares with similar properties in both long and 
short bursts, as well as accounting for the observed corre- 
lation between the flare duration and its arrival time (see 
also Piro & Pfahl 2006). Other disk-based models involve 
magnetic instabilities (Proga & Zhang 2006; Giannios 2006). 
Staff, Ouyed & Bagchi (2006) proposed instead a mechanism 
based on state transitions in the quark-nova scenario, while 
the role of a magnetar was discussed by Gao & Fan (2006) 
and Cea (2006). 

4.1 Refreshed shocks 

If an "internal" flare, i.e., a flare due to internal dissipa- 
tion, is observed in a GRB, it must be followed at later 
time by an external flare. The material responsible for the 
internal flare keeps expanding at constant speed, and must 
eventually catch up with the external shock, which deceler- 
ates continuously. It may however be difflcult to observe the 
second flare, since it may be weak and dispersed over sev- 
eral orders of magnitudes in time. If, however, the second 
flare can be singled out and flrmly associated to an external 
flare, it would provide additional important information for 
our understanding of the dynamics of the system. 
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Consider a flare system made by an internal flare begin- 
ning at time tint followed by a second, broader, flare starting 
at time text, where text S> tint. The external flare time can 
tell us the Lorentz factor of the material producing the in- 
ternal flare: 

where E is the isotropic equivalent kinetic energy of the flre- 
ball integrated up to the time of ejection of the late shell, 
n the numeric density of the external medium, and nip the 
proton mass. We have assumed a uniform external medium. 
Equations can be easily generalized to an arbitrary density 
profile, but will be simple only in the case of adiabatic evo- 
lution. Equation (|10|l depends only on measurable quanti- 
ties. An example of a solvable system is the flare system of 
GRB 050502B. Falcone et al. (2006) were able to identify 
an initial strong flare and a possible refreshed shock flare 
at later times. They were able to measure the bulk Lorentz 
factor of the late engine ejecta to be Foj ^ 20, ejected from 
the inner engine approximately 300 second after the onset. 
By itself this observation only tells us that the late ejections 
are likely to have a smaller Lorentz factor than the early 
ones. Solving a large numbers of flare systems in different 
GRBs would however provide us with an ensemble of late 
ejections that would allow us to constrain both the proper- 
ties of the engine and of the dissipation/radiation mecha- 
nism. The knowledge of the bulk Lorentz factor of the late 
shells would provide us with a sample of flares for which 
the comoving spectra are known, ridding the prompt emis- 
sion models of one of the unknown parameters: the Lorentz 
factor. 



5 SUMMARY 

We have analyzed the timescale of flares produced by several 
mechanisms both due to internal dissipation and external 
shock phenomena. We have shown that all prominent exter- 
nal shock flares occurring before the jet break must have a 
long time scale 5t ^ t. We also showed that internal dissipa- 
tion taking place in slow shells ejected immediately after the 
end of the prompt phase must satisfy a similar constraint: 
5t ^ 0.25t. The only mechanism capable of producing flares 
with a short time scale as observed by Swift is identified with 
late activity of the inner engine. As a consequence, the de- 
tection of X-ray flares at times as long as 10^ seconds after 
the GRB onset implies that the inner engine is active for at 
least 10^ s. This has important implications for the fueling 
of the engine and the ejection mechanism. 
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